Signal to Noise Ratio Estimation of the ASCENDS 
CarbonHawk Experiment Simulator (ACES) for 
Atmospheric C0 2 Measurement 


Songsheng Chen 1 , Larry B. Petway 1 , Yingxin Bai 2 , Fenton W. Harrison 1 , 
Joel F. Campbell 1 , Byron L. Meadows 1 , Bing Lin 1 , 

Michael D. Obland 1 , and Edward V. Browell 2 


J NASA Langley Research Center, 5 N Dry den Street, Hampton, VA 
23681 

2 Science Systems & Applications Inc.,1 Enterprise Parkway, 
Hampton, VA 23666 

3 STARSS II Affiliate NASA Langley Research Center 


SPIE REMOTE SENSING 

Toulouse France, Sept. 21-24, 2015 


Outlines 


> Introduction 

> Maximum Likelihood Estimation (MLE) for Multiple 
Linear Swept-Frequency Sine- Wave Detection 

> An Airborne Multiple Swept-F requency (SF) Intensity 
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• MLE of multiple swept-frequency digitized IM-CW 
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C0 2 Sources and Sinks (Missing?) 

Fossil Fuel Emissions of C0 2 and Atmospheric Buildup, 1958-2008 
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Maximum Likelihood Estimation (MLE) -1 



> MLE for analysis of the parameters of multiple-frequency 
digitized sine-wave signals based on a statistical model, 
all the digitized signals, sin) , is sum of the expectation 
value of the signals and Gaussian-distributed random 
noise with the probability, 
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E(s{n)) is expected value of s(n) and is variance of the 
noise. 

> MLE formula 
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Maximum Likelihood Estimation (MLE) -2 


> Maximum Log-Likelihood Estimation Function / 
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Maximum Likelihood Estimation (MLE) -3 


> Multiple Digitized Sine Wave Signals (numbers of k) 
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Maximum Likelihood Estimation (MLE) -4 


> Maximum Log-Likelihood Estimation Function l for 
Multiple Digitized Sine Wave Signals with constant 
Frequencies 
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This is the formula for the computer simulation of the 
MLE of the airborne multiple swept-frequency IM-CW 
ACES system 
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An Airborne Multiple SF IM-CW ACES 


> Advances of airborne technologies in support of NASA 
ASCENDS mission 


Airborne Platforms: 



o High and large altitude range: ~1.0 km - ~20 km 


EDFA-based high-power Laser 
transmitter: 


up to 3x1 OW 
(Avg.) 
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HgCdTe APD Array-based high-sensitive 
detector: 

NEP: 2.4fW/Hz 1/2 
Temperature < 77k 
Bandwidth: 

~4.9 MHz @ gain 10E( 
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An Airborne Multiple SF IM-CW ACES 


> Airborne multiple SF IM-CW ACES operated at Three 

Different Wavelengths around C0 2 Absorption Line of 1571nm 

Two-Way Optical Depth nline Weighting Functions 

at C0 2 Absorption Line at of 1571nm 



Calculated based on HITRAN 2008 Database, US standard atmosphere 

>-om => -<f 5(, P m ; '•off2= > -o+ 50 P m 
>^,=1571. lllnm —Center Wavelength 
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An Airborne Multiple SF IM-CW ACES 

> EDFA-Based high-power laser transmitter 
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An Airborne Multiple SF IM-CW ACES 


> Multiple SF IM-CW ACES optical and electrical digitized 
receiver 
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An Airborne Multiple SF IM-CW ACES 


> LIDAR Equation for the airborne multiple SF IM-CW ACES a 
three different aavelengths 

• Total Received Signals under the Intrinsic Limit of Shot Noise 
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(*) 3,° ^ • Laser signals, Signal amplitude, Modulation index at Xi ; 
z\ and OD a . Two-way atmospheric transmittance and optical depth at Xi ; 

7j q and rj opt : Detector quantum efficiency and system optical efficiency; 

D t , r s , r a eff , and (A X) : Dia. Of Telescope, Surface reflectance, Effective 
atmospheric transmittance and bandwidth near C0 2 absorption line at 1571 nm 
R, J, 0, (AX): Transmitter-Target range, Solar background induced photon 
radiance and scattering, and FOV of telescope 
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An Airborne Multiple SF IM-CW ACES 


> Maximum Likelihood Estimation of the Digitized Multiple- 
Frequency IM-CW Signals at Three Different Wavelengths 

• Gaussian noise assumption: Detected laser signal photons, 
solar background photons, and detector dark current - 
contributed numbers of electric charges (electron) during 
one sampling period are more than 15-20 

Total detected signal in a simple format: 

S(r) = S° c + X Sf cos(2 7f t m t + <p)+ N g (?) 

1=1 

= .S’ ( ° + X (s; cos(2 nfrt)+S; sm(2nfrt))+N G {t) 

i = 1 
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An Airborne Multiple SF IM-CW ACES 



> Maximum Likelihood Estimation of the digitized multiple SF 
IM-CW Signals at three different wavelengths 

• Gaussian noise assumption: Detected laser signal photons, 
solar background photons, and detector dark current - 
contributed numbers of electric charges (electron) during 
one sampling period are more than 15-20 

• Gaussian noise with constant variance: Simpler case for 
solar background photon and detector dark current 
contributed shot noise is much larger than signal shot noise 

cr 2 = a 2 proportional to total detected solar background 
photons at low detector dark current. Simple MLE formula 
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An Airborne Multiple SF IM-CW ACES 


> Maximum Likelihood Estimation of the digitized multiple SF 
IM-CW Signals at three different wavelengths 
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> Simplified linear equations 

Coefficient matrix R: off-diagonal components can be zero under 
following cases, 

* Sampling frequency is much larger than modulation frequencies 

** Specified sampling frequency, all modulation frequencies, and 
integration time periods. 

Our simulation condition: 

SF signals at three different wavelengths: 

Sampling frequency: 2 MHz 

Swept frequencies: 100.4 kHz-600.4 kHz; 105.2 kHz-605.2 kHz; 
and 110.8 kHz-610.8 kHz 

Frequency increase rate (chirp rate): 2500 MHz 

Orthogonal sweeping time: 5 ms (25 swept ) 

Integration time periods: 5ms, 25ms, 50ms, 100ms 
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An Airborne Multiple SF IM-CW ACES 


> The Fractional/Relative Random Error of the C0 2 measurement 
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Fractional/relative errors of the C0 2 column density measurement 

Differential optical depth of the transmitted signals in the atmosphere 
between on-line and offline (1 or 2) wavelengths 

Signal to Noise Ratio (SNR) at on-line and off -line (1 or 2) wavelengths 
on the transmitter detector for laser power monitor 

Signal to Noise Ratio (SNR) at on-line and off-line (1 or 2) wavelengths 
on the detector of the receiver for signal measurement 

r o„ and^ are two-way transmittance at on-line and off-line (1 or 2) wavelengths 
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An Airborne Multiple SF IM-CW ACES 


System Parameters for Computer Simulation of MLE 

Altitude: 2-15 km 

Receiving telescope: 

Effective Diameter: 3*(0.17/V3)m 

(3 7” Ritchey-Chretien telescopes with center blocks in ~2” dia.) 



FOV (full angle): 

495 prad 

Over all optical transmittance 

8.5% 

Optical filter bandwidth: 

2.7 nm 

Laser transmitter: 

Beam divergence: 

300 prad 

Modulation index: 

0.9 

Photodetector: (HgTeCd APD 8*8 array) 

Noise Equivalent Power (NEP): 

2.4fW/Hz 1/2 @~77K 

Quantum efficiency: 

0.8 

Ground radiance near 1571 nm: 
(W/m2/sr/jim) 

ocean/vegetable 

1. 7/5.0 

Surface reflectance: 

ocean/vegetable 

0.08/0.31 

Atmospheric transmittance for radiance 

95% 
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Relative error of C0 2 column density 


An Airborne Multiple SF IM-CW ACES 

Computer Simulated MLE Results 
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Total Power distributed to different wavelength based on the atmospheric transmissions 

9/13/2015 


Document Title 


22 


SNR at on-line wavelengths 



Conclusions and Discussions 



> MLE of multiple digitized swept-frequency (SF) 
Intensity wave signals with Gaussian-distributed 
noise 

> SNRs and relative errors of the C0 2 column density 
measurement from airborne at different altitudes 
have been simulated for the airborne SF IM-CW 
ACES 

> SNRs increases as integration time increases up to 
0.1s (linearly in log scales) 
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